Generation and detection of a spin entanglement in nonequilibrium quantum dots 
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Spin entanglement between two spatially separated electrons can be generated in nonequilibrium 
interacting quantum dots, coherently coupled to a common lead. In this system entangled two- 
electron states develop which are Werner states with an imbalance between singlet and triplet prob- 
abilities. We propose a multi-terminal, multiply-connected setup for the generation and detection 
of this imbalance. In particular, we identify a regime in which the formation of spin entanglement 
leads to a cancellation of Aharonov-Bohm oscillations. 
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I. INTRODUCTION 

The controlled generation and detection of entangle- 
ment of quantum states remains one of the fundamental 
challenges of quantum physics. Experiments with en- 
tangled photons have already entered the realm of ad- 
vanced quantum communication and cryptography. In 
solid-state systems electron spins are considered as prime 
candidates for the demonstration of such effects. Spin de- 
grees of freedom are only weakly coupled to the environ- 
ment, which leads to long decoherence times [l], and 
coherence lengths well exceeding the micrometer range. 
Several setups involving quantum dots have been pro- 
posed to create and detect pairs of spatially separated, 
spin-entangled electrons. The schemes rely, e.g., on the 
extraction of Cooper pairs from a superconductor Q , or 
the separation of a pair of entangled electrons from a sin- 
glet ground state of single 0, Q or double quantum dots 
jg, 0] • Evidence of entanglement in these systems can be 
obtained from noise measurements Q or coincidence de- 
tection 0. In Ref. @ a detection scheme was proposed 
where pure singlet and triplet states can be distinguished 
by Aharonov-Bohm interferometry. 

In this article we analyze a setup with two spatially 
separated quantum dots, illustrated in Fig. [T] which al- 
lows the creation and detection of entanglement of spa- 
tially separated electron spins. The two dots are cou- 
pled coherently to a joint source electrode on the left and 
to two independent drain reservoirs at the top and bot- 
tom. When a bias voltage is applied, electrons are driven 
from the source via the dots to the drain electrodes. A 
nonequilibrium, mixed electronic state is created. For 
appropriate values of the voltage, due to the strong on- 
site Coulomb repulsion, a state with two electrons (one 
on each dot) has a high probability. It turns out to be of 
the form of a Werner state, with a strong enhancement of 
the singlet component, although singlet and triplet states 
are degenerate [Icj . 

In order to detect the entanglement, an additional joint 
reservoir, closing the Aharonov-Bohm geometry, is cou- 
pled to the system. The current to this reservoir is stud- 
ied in the cotunneling regime. Under certain conditions 




FIG. 1: Two quantum dots (u and d) coupled to a joint 
reservoir on the left hand side (L) and to two separate leads at 
the top and the bottom (U and D). The second joint reservoir 
on the right (R) closes the Aharonov-Bohm ring and is used 
to probe the state of the system. 



the entanglement leads to a suppression of the Aharonov- 
Bohm oscillations. Hence this part of the setup serves as 
probe of the state of the system and as detector for en- 
tanglement. Exploring the results in a wide parameter 
range we identify regimes where the double dot has a 
large probability to be in a singlet state. The predicted 
behavior provides a proof of concept for the entangle- 
ment generation in coherently coupled, nonequilibrium 
quantum dots. 



II. MODEL 

We consider the system shown in Figure [U It consists 
of two single- level quantum dots, i = u, d, with strong on- 
site Coulomb repulsion U described by the Hamiltonian 
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The dots are coupled to joint reservoirs on the left and 
right forming an Aharonov-Bohm loop, and to two sep- 
arate leads "up" and "down". The total system is then 
modelled by H = i/dots + H Tes + H t with the four leads, 
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r = L, R, U, D, with spin degenerate energies e r k, 

Hrcs — ^ &rk ^ rkcT ^rk(T , (2) 
rka 

and the four tunnel Hamiltonians H t — H t L+Htn+H t Tj + 
H t u- Tunneling between the upper dot and reservoir is 
simply described by 

H t V = ^2 [tu A,a a Uka + h.C.] , (3) 
ka 

and similar for the lower pair. Tunneling between the 
joint leads on the left and the right, on the other hand, is 
sensitive to the phases of the matrix elements. Without 
loss of generality we can distribute them in a symmetric 
gauge, 

H th = [(*L U e- iv/4 4, + i Ld e +W4 4 CT ) a Lka + h.c] , 

ka 

H t R = [(*Rue +i ^ /4 cL + t Rd e~^/ 4 cl^ a Rka + h.c. 

ka 

(4) 

All tunneling matrix elements arc taken to be indepen- 
dent of spin and energy, and they can be chosen real and 
positive. The leads are assumed to be in local equilib- 
rium, described by electrochemical potentials [i r . 

To keep the discussion transparent, we assume some 
symmetries. We set the upper and the lower reservoir to 
the same potential, /iu/D — Mu = Md, and restrict our- 
selves to degenerate dot levels, e = e u = e^. We, further- 
more, assume symmetric coupling strengths of the upper 
and lower dot to the left reservoir, t^u = ^Ld = th, and 
similar for the right reservoir, t Ru = i^d = £r, and we 
assume % = to = £u/d- Defining the tunneling strength 
by IV = 27rt^ N r , where N r is the density of states of the 
relevant reservoir, we are left with three parameters IY, 
Tr, and T v / D , which gives us enough freedom to tune the 
system into the interesting regimes. Double occupancy 
of each dot is suppressed by a strong Coulomb repulsion 
U, for simplicity assumed here to be larger than all other 
energy scales. However, the generalization to finite U is 
straightforward. We apply two bias voltages measured 
relative to /il = 0: a strong bias voltage between the left 
and the upper and lower reservoirs, eVu/D = Mu/D: and a 
small transport voltage between the left and right reser- 
voir, eVa = ^r, used to probe the state of the system. 

In essence, the preparation and measurement scheme 
analyzed in the following is based on two observations. 
First, the components of the density matrix which de- 
scribe the doubly occupied system with one electron on 
each dot is a so-called Werner state, a mixture of a sin- 
glet and three equivalent triplet states. The distribution 
between singlet and triplet probabilities can be tuned in 
a certain range by the coupling strengths and the poten- 
tials applied to the system. For illustration we consider 
the state obtained immediately after charging the double 



dot with two electrons with opposite spin from the left 
reservoir 

(e-^cl + e+^cl)(e-^ci s + e +i ^ 4 4)|0) 
= e-^ 2 \aa, 0) + e +iv / 2 |0, aa) + \a, a) - \a, a) . 

The minus sign in the last term is a direct consequence 
of the Fermi statistics. In view of this sign and the Pauli 
principle it is obvious that the equivalent expression for a 
pair of electrons with equal spin vanishes. A strong onsite 
Coulomb interaction prohibits double occupancy of each 
dot. Hence the state reduces to a perfectly entangled, 
pure singlet state 

\S) = (\a,a)-\a,a))/V2. 

In equilibrium, however, the state has relaxed to a 
completely unentangled mixed state with uniformly dis- 
tributed probabilities for singlet and triplets. 

The second observation is that for the cotunncling 
rates the phases of the transport contributions for sin- 
glet and triplet states are shifted relative to each other 
by 7r. Hence they yield different interference patterns in 
the Aharonov-Bohm interferometer. 

The strategy, therefore, is to use the subsystem com- 
posed of the left, upper, and lower reservoirs, and the 
dots to prepare the state of the system. In particular, 
we tune the ratio of singlet and triplet probabilities via 
the bias voltage V\j /d ■ The signatures of this state in the 
transport are monitored by the Aharonov-Bohm interfer- 
ometer, which is operated in the cotunneling regime for 
the tunneling to the right reservoir. The main message 
of this article is that for a setup with cascading coupling- 
strengths, Tl ^> Tu/pj ^ Tr , an enhanced probabil- 
ity for singlet over triplet formation is associated with a 
suppression of the oscillations in the conductance of the 
Aharonov-Bohm probe. 

III. PREPARATION OF STATIONARY 
ENTANGLED STATES 

We investigate the system in the frame of the real-time 
diagrammatic technique Refs. [ll|, [H, [l3| which allows us 
to go beyond an orthodox master equation approach. In 
particular we account for the quantum coherent time evo- 
lution of the double dot and focus on correlations between 
the electrons. The elements of the reduced density ma- 
trix are given by = (x'|tr rcs p|x), where the reservoir 
degrees of freedom of the full density matrix are traced 
out, and x and x' label general many-body states of the 
double dot. For strong Coulomb repulsion, the reduced 
Hilbert space is spanned by nine basis states |x UJ Xd) with 

Xi = o,u. 

To describe the symmetries implied by tunneling we in- 
troduce a tailored basis. The state of the empty system 
and its probability are denoted by |0,0) and po, respec- 
tively. A single electron of spin a can occupy, with overall 
probability for single occupancy p\ — ^2 irT Pia , either the 



3 



upper or the lower dot, or, in general, any coherent super- 
position of these two. Thus, the corresponding block of 
the density matrix describes a quantum two-state system, 
which we express an isospin for each spin component, I a . 
A natural basis are the two states with one electron in 
the upper or in the lower dot. However, in this basis 
the joint left and right reservoirs couple via the Hamil- 
tonian ((H) to the non-diagonal components I a ■ hil/r, of 
the isospin, with 

m L /R = (cos - , ± sin - , 0) , (5) 

which introduces a non-diagonal time evolution of the re- 
duced density matrix. Electrons entering the double dot 
from the left reservoir generate an isospin polarization 
in direction of niL, whereas electrons entering from the 
right are polarized in the direction of mR. We further 
note, that due to the spin symmetry of the Hamiltonian 
the two spin components are completely equivalent. We 
will formulate equations for the expectation values of the 
isospin for either spin component, (If) = (I|) = 
To simplify expresssions we dropped in the last form the 
angular brackets. 

The doubly occupied subspace with one electron on 
each dot is naturally spanned by the spin singlet |S), 



which we find with probability ps , and the three triplet 
states |T + ), |T ), and |T_). Due to the spin symmetry 
of the Hamiltonian the three triplet states are equivalent 
and are occupied with equal probability pt — Pt± — 
Pt/3. The corresponding block of the density matrix 
can be represented in the Werner form [3] 

W(F) = F |S){S| + (1 - F) l4 "] S)(S| , (6) 

where the coefficient F = ps/(ps + Pt) denotes the 
Werner fidelity. Values in the range of 1/2 < F < 1 
indicate a quantum correlated mixed state from which 
arbitrary entanglement can be distilled by suitable pu- 
rification protocols [HI, [l6| • As we will show below, the 
Werner fidelity of the stationary state can be tuned in a 
certain range by the tunnel couplings and the applied bias 
voltage. Regimes with F > 1/2 are feasible if the double 
dot is predominantly charged from a common reservoir. 

The kinetic equations for the reduced density matrix 
can be arranged in vector form with the diagonal prob- 
abilities collected in p = (j>o,Pi,PStPt) and the state 
of the singly occupied subspace specified by the isospin 
lyi = (I x ,I y ,I z ). In lowest order in the coupling 
strengths Tl^r, Tu/d we obtain 
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(7) 



r=L,R 



/>..+ (/r-i)f + (l-/r)f -(l-/r)| 



m r - ^ r r (1 + fr) It/1 

r=L,R,U/D 



where the Fermi distributions f r = 1/ (e /3 ^ e_ + 1) in the 
respective reservoirs are evaluated at the dots' energy e. 
The Equations ([7]) extend our earlier work (To| , where we 
considered a double-dot system coupled to three reser- 
voirs, but without the Aharonov-Bohm probe and right 
reservoir. 

To generate and detect a singlet-triplet imbalance, the 
system parameters have to be chosen specifically. First, 
the tunnel coupling to the right lead should be much 
weaker than the other ones, Tr <C Tl^u/Dj hi order to 
ensure that the state of the double dot is not affected by 
the measurement via the interferometer. Furthermore, 
only a small bias voltage is applied across the Aharonov- 
Bohm ring, and the dot levels are kept off resonance such 
that the flux-sensitive linear conductance is dominated 
by cotunneling processes. Finally, in order to find with 
large probability two electrons in the system, the rate 
for charging should be much larger than for discharging, 



r L » r u/D . 

In the following, we compare two situations: 

(i) For vanishing bias voltage between the left and up- 
per/lower reservoirs, the double dot remains in equilib- 
rium, and F = 1/4. This defines the reference for com- 
parison. 

(ii) By applying a high bias voltage between the left and 
upper/lower reservoirs, we generate a singlet-triplet im- 
balance with F > 1/2. 



IV. PROBING THE SINGLET-TRIPLET 
IMBALANCE 

Cases (i) and (ii) lead to significantly different interfer- 
ence signatures in the transport through the Aharonov- 
Bohm ring. While the first one shows a strong flux de- 
pendence and large amplitude of the Aharonov-Bohm os- 
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dilations, they are suppressed in the second case, where 
the Werner fidelity approaches 1/2 if the double dot is 
predominantly occupied with two electrons. To demon- 
strate this difference we calculate the linear conductance 
through the Aharonov-Bohm interferometer. The contri- 
bution to first order in the coupling strength is evaluated 
within the diagrammatic technique. It turns out to be 
small due to our choice of parameters (compare also Fig- 
ure [3]) . For the dominant second-order contribution we 
calculate the stationary density matrix according to ([7]) 
and obtain the cotunneling rates from Fermi's Golden 
rule pi 13113 



(2ir) 2 N-^Ni 



a 2^ 



(Xf k { \H t \n)(n\H t \xi k 



71 



LO - E„ 



(8) 



Here the initial double-dot state is denoted by xu an d all 
possible final states are summed over. The energy E n of 
the intermediate virtual state \n) is either 0, e or 2s, and 
lu is the energy of the incoming lead electron. 

Gathering all cotunneling terms corresponding to the 
four scenarios where the dot system is empty (% = 0), oc- 
cupied with a single electron with isospin Ij n (x = 1); or 
occupied with two electrons either in a singlet (x = S) or 
a triplet (x = T) state, we find the second order contri- 
bution to the linear conductance Gr = dl^/dV^y _ i 



2 f 

=~g{ V )Ke J 



TtT 



Ll R 



■/»• (9) 



The Aharonov-Bohm flux dependence is included in the 
dimensionless conductance g(ip) = X) x 5x(^)' w ith 

5o M = 2(1 + cos tp)p ((f) (10a) 
gi((p) = 2 (3 - cos if) pi(<p) + 16 (m L - m R ) • I T/i 

(10b) 



.'/Slr"'l ( 2 ~ \ C0Sl P j > 



1 1 



cos</3 pr(y) 



(10c) 
(lOd) 



The flux sensitivity of the conductance is twofold. First, 
the phase factors from the tunneling Hamiltonian are 
taken into account in the coherent summation of pro- 
cesses in the cotunneling rates (JSJ). In particular, as 
discussed in [9j, the phase dependence of the contribu- 
tions from an initial singlet and triplet state (Equations 
(|10cj) and (JTOdJ)) are shifted relative to each other by 
7r. We emphasize that these results rely on the degen- 
eracy of singlet and triplets. All four constitute equiva- 
lently possible final states of cotunneling processes with 
an initially doubly occupied system. If there was a 
large singlet-triplet energy splitting with transitions be- 
tween singlet and triplet states being suppressed, then 
5s,t(v) = (1 + cos (^)ps.t/4 without any distinction be- 
tween the contributions. 

The second kind of flux dependence results from the 
stationary state distribution of the system. In general 



nonequilibrium situations, the reduced density matrix, 
represented by po, pi, Psi an d Pt, is influenced 

by the Aharonov-Bohm probe and becomes phase sensi- 
tive. This effect is weak for weak tunnel coupling Tr <C 
Tl, Tu / D where the state of the system is only weakly in- 
fluenced by the AB-probe. In the limit Tl ^> Tu/d 3* Tr 
and the double dot charged from the left reservoir the 
isospin points approximately in direction of niL such that 
gx(<p) w 2(3- cos<^)pi - 16 (1 - cosf) |I T/i |. 

We compare now the two scenarios (i) and (ii). For 
low dot energies in an equilibrium situation, s <C Hr, the 
double dot is mainly occupied with two electrons with 
p T w 3/4 and ps ~ 1/4. As a consequence, in this equi- 
librium reference situation there is a good visibility of the 
Aharonov-Bohm oscillations, 



5 (i) M 



1 1 



■ COS if . 



(11) 



In contrast, for strong bias voltage, the phase dependence 
is suppressed. To see this, we expand the density matrix 
for a strong asymmetry in the coupling to the left, and 
the upper and lower reservoirs, x = T^j/d/IY <C 1, and 



0(x 2 ), pi = 4a;/3 + 0(x 2 ), I T/i w x/6 m L + 
= 1/2 - x/2 + C(a; 2 ), and p T = 1/2 - 5x/6 + 
0(x A ). Thus, the dimensionless conductance approaches 



obtain po - 
0(x 2 ), p s 



14 29 
y-^cos^ 



(12) 



For strong asymmetry, Tu/d *C Tl, the oscillations of the 
singlet and triplet terms cancel each other and the phase 
dependence of the conductance vanishes. This suppres- 
sion of the Aharonov-Bohm amplitude allows us to detect 
the singlet-triplet asymmetry generated in system. 

We remark that if the system is prepared in a state 
with maximal entanglement, p$ = 1 or px — 1 5 as consid- 
ered in Ref. 0] , then the visibility of the Aharonov-Bohm 
oscillations would be even twice as large as in equilibrium, 
g = | =F \ cos if. For a steady-state scenario considered 
in this paper, however, it is not possible to have both a 
large overall probability of double occupancy of the dou- 
ble dot and a maximal singlet-triplet imbalance at the 
same time. To optimize the probability for double occu- 
pancy we choose Tl ^> for which p$ — pt = 1/2, 
which indicates finite but not maximal entanglement. 

The plots in Figure [5] display the linear conductance up 
to second order in the coupling strength based on a full 
solution of the stationary equations . If the system is 
close to equilibrium, |eVu/D| ^ k^T, the triplet probabil- 
ity is large (F w 1/4) , and the conductance is dominated 
by a positive cos tp kind of oscillation (compare Equa- 
tion (fTT| ). In contrast, the interference is suppressed if 
the system is driven into a singlet-triplet imbalance by 
charging the double dot from the joint left lead, F ~ 1/2. 

A finite singlet-triplet relaxation reduces the imbal- 
ance between singlet and triplet in the stationary state, 
and the Aharonov-Bohm oscillations are restored. To 
estimate the influence of spin flip and dephasing we in- 
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bias voltage, however, its phase dependence is weak since 
for strongly asymmetric coupling the isospin is parallel to 
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FIG. 2: The linear conductance of the Aharonov-Bohm sub- 
system up to second order in the coupling strength is plot- 
ted versus the Aharonov-Bohm phase. The parameters are 
r L = IknT, r u/D = 0.1 kuT, r R = 0.01 fc B T, e = -10fc B T 
and |eVEt| <C fc B T. For strong bias voltage between the 
left and the upper/lower reservoirs, eVu/D = — 20fc B T (left 
plot), and in equilibrium (right plot), the ideal system (i) 
is compared to results for singlet-triplet decay rates of (ii) 
Tst = 0.5r u/D and (in) Fst = 10r u/D . 
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FIG. 3: The individual contributions to the linear conduc- 
tance are plotted with and without singlet-triplet relaxation. 



The first order, and second order terms, Gr' x , are spec- 

ified for equilibrium (right plots) and for strong bias volt- 
age, eVu/D = — 20feT (left plots). The parameters are 



-,(2) 



i l = 1 k B T, r u/D 
and leVy < k B T. 



o.i/cbT, r 



R 



0.01 k B T, e 



-10fc B T 



troduce phenomenological transition rates between sin- 
glet and triplets. For simplicity we choose all of them 
equal to Tst- To observe a significant suppression of 
the Aharonov-Bohm oscillations a relaxation rate smaller 
than the current rate between source and drain is re- 
quired, i.e. Tst < Tu/d m our case. 

In Figure[3]the first and second order conductance con- 
tributions are plotted individually. The transport sup- 
ported by the equilibrium state is dominated by the sec- 
ond order triplet contribution and the Aharonov-Bohm 
oscillations are well visible. The singlet-triplet relaxation 
does not affect the behavior. If the electrons are strongly 
driven into the double dot from the joint left reservoir, 
a singlet-triplet imbalance forms and the oscillations of 
the singlet and triplet contributions cancel each other. 
A finite relaxation rate partially destroys the imbalance 
and the cancellation. The conductance supported by the 
singly occupied states becomes quite strong for a large 



CONCLUSIONS 



Here we proposed a setup for the generation and de- 
tection of spin entanglement between two spatially sepa- 
rated electrons. A bias voltage applied across a double- 
dot system with strong onsite Coulomb repulsion, co- 
herently coupled to a joint source electrode and individ- 
ually coupled to separate drains, creates an imbalance 
between singlet and triplet probabilities. An asymmetry 
in the coupling of source and drain increases the overall 
probability to find two excess electrons in the double dot 
system and for them to form a Werner state. The un- 
derlying mechanism of entanglement generation relies on 
the specific nonequilibrium situation and is fundamen- 
tally different from those schemes in which two electrons 
are extracted from a singlet ground state. 

To detect the singlet-triplet imbalance we added an 
Aharonov-Bohm probe, obtained by coupling a second 
joint reservoir coherently to the double dot. The enclosed 
flux controls the interference pattern which depends on 
the state of the double dot. In particular, the cotunneling 
transport through the Aharonov-Bohm ring is sensitive 
to an imbalance in the singlet-triplet distribution. 

We compared two situations: For vanishing bias volt- 
age between the joint source and the drain reservoirs the 
double dot is in equilibrium, and singlet and triplet are 
uniformly distributed, and the Aharonov-Bohm conduc- 
tance is strongly flux dependent. On the other hand, 
by applying a strong bias voltage a singlet-triplet imbal- 
ance is generated, and the Aharonov-Bohm oscillations 
are suppressed. By choosing asymmetric couplings we 
can ensure that the double dot was predominantly occu- 
pied with two electrons. 

For an experimental realization one needs a geome- 
try with spatial separation of the dots shorter than the 
phase coherence length. Finite decoherence due to spin- 
orbit effects, hyperfine coupling or coupling to some other 
bath will reduce the overall visibility of the Aharonov- 
Bohm oscillations. Nevertheless, the extra suppression of 
the Aharonov-Bohm oscillations for strong bias voltage 
as compared to equilibrium will indicate a singlet-triplet 
asymmetry. The probability to find the system occu- 
pied with two electrons can be enhanced by tuning the 
dot levels well below the Fermi energy of the probe and 
choosing a strong asymmetry in the coupling of source 
and drain. For the detection of a suppressed phase de- 
pendence caused by a singlet-triplet imbalance a tunnel- 
ing rate between source and drain much larger than the 
spin decoherence time is required. 
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